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We propose an electronic model for the recently discovered hole doped compound Y2- x Ca x BaNiOs. 
From a multiband Hamiltonian with oxygen and nickel orbitals, a one band model is derived. Holes 
are described using Zhang-Rice-like S=l/2 states at the nickels propagating on a S=l spin chain. 
Using numerical techniques to calculate the dynamical spin structure factor S(q, uf) in a realistic 
regime of couplings, spectral weight in the Haldane gap is observed in agreement with neutron 
scattering data. Low energy states with S=3/2 appear in the model. Several predictions are made 
to test these ideas. 
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Spin-liquid ground states in Heisenberg S=l chains 
and S=l/2 ladders have been studied theoretically as 
paradigms of disordered nonclassical systems. |l],|2| These 
spin models can be physically realized in several com- 
pounds. An important issue is the effect of hole doping 
on these systems. Theoretical studies of doped S=l/2 
ladders have shown that the spin gap survives in the pres- 
ence of holes, and that the spin-gapped phase is favorable 
for superconductivity. Behavior indicating a spin gap 
has also been observed experimentally in some under- 
doped high-Tc cuprates. However, the case of the doped 
S=l chains has been considered only very recently in the 
context of the S=l metal oxide Y2BaNiOs. || Lightly 
doping this compound with Ca introduces hole carriers in 
the chains. Two remarkable experimental features were 
observed upon doping, namely (i) the reduction of the re- 
sistivity pdc by several orders of magnitude leading to a 
one dimensional conductor, and (ii) the creation of states 
inside the Haldane gap as revealed by inelastic neutron 
scattering (INS) data. {§ 

To understand these in-gap states, spin systems with 
site or bond impurities have been recently proposed. ||] 
Holes are assumed to be so strongly localized that their 
mobility is neglected. In another approach dynamics was 
provided to the holes || but the Ni-Ni exchange J was 
assumed larger than both the hole hopping amplitude 
and the short-bond NiO exchange. While studying these 
limits is instructive since gap states are produced, the 
drastic reduction in pd c observed experimentally upon 
doping Q suggests that holes are likely mobile over sev- 
eral lattice spacings. Q Thus, in this paper we propose a 
new model for Y2- x Ca x BaNi05 with fully mobile S=l/2 
holes interacting with S=l spins, which is studied using 
realistic values for the couplings. 

In Ni ++ surrounded by oxygens, d 3z 2_ r 2 and d x 2_ y 2 
are the two active orbitals if deviations from the perfect 
octahedral symmetry in NiOg are considered. || Then, as 
a Hamiltonian for the Ni — O chains in the hole notation 
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i (j) denotes Ni (O) sites, a — 1(2) corresponds to the 
Ni orbitals d 3z 2_ r 2 (d x 2_ y 2). The Coulomb repulsion is 
Ud (Up) at the Ni (O) sites, and A is the charge-transfer 
energy. di CTCt are hole operators corresponding to a Ni 
site, spin a and orbital a, while pjo- are oxygen hole 
operators. The last term in Eq.(l) is a ferromagnetic 
coupling between the Ni holes on different orbitals (us- 
ing Si Q = d; Q (Tdi Q /2), which enforces Hund's rule. This 
term is important to produce the expected S=l state 
in Ni ++ . From an analysis of various spectroscopies for 
Li x Nii_ x O, it was found that A = 6.0eV, U d = 9.5eV, 
and Up = 4.6eV. M Since Ud > A then the compound 
is in the charge-transfer regime. ]io| Jfjund is obtained 



from the energy levels of a (NiOg) lu ~ cluster. |llj The 
difference in energy (|JHund|) between the S=l and 
the S=0 1 E levels is 1.3eV. The hopping amplitudes are 
tpdi = 1.3eV and t pd 2 = 0.75, according to a cluster cal- 
culation. 0^3) 

The study of mobile holes in model Eq.(l) is a difficult 
task since the energy scale is eV, while the interesting 
physics for Y2BaNiOs occurs at the J-scale of about 30 
meV. Thus, we need to construct an effective low-energy 
Hamiltonian from Eq.(l) which here we carry out fol- 
lowing the work of Zhang and Rice in their reduction of 
the 2D multiband Hubbard model to the t-J model. JlJ] 
In Y2- x Ca x BaNi05 holes populate oxygens, |3| and they 
have a strong exchange J' with neighboring Ni. This may 
lead to low energy states with S=l/2 or S=3/2, after the 
oxygen S=l/2 hole mixes with neighboring S=l spins. 
To study this effect, we solved a small cluster O-Ni-0 de- 
scribed by Eq.(l) in the three holes subspace (two of them 
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producing the Ni S=l, and the other providing an extra 
oxygen hole). From the spectrum, and using the pa- 
rameters given before, we indeed found a S=l/2 ground 
state, with the S=3/2 state located ~ 1.3 eV higher. |l5fj 
This energy difference is robust even if the parameters 
of Eq.(l) are modified within acceptable windows to ac- 
count for experimental uncertainties. Then, proceeding 
a la Zhang-Rice we project out all the states of the three 
holes O-Ni-0 cluster but the Ni-centered S=l/2 ground 
state. [|l6| Eq.(l) thus reduced to an effective model with 
states only at the Ni sites which can be S—l ("spins") 
or S=l/2 ("holes"). Thus, in our model the study of 
carriers on Y2_ x Ca x BaNi05 amounts to the analysis of 
S=l/2 hole-like states in a S=l background. ]l7| 

To derive the allowed hole hopping processes it is 
important to remember the composite character of the 
S=l/2 hole state. To guide the intuition a graphical rep- 
resentation is useful (Fig. la). The low energy states with 
S=l and 1/2 located on nearest neighbors Ni sites are ac- 
tually represented as five holes (four in the Ni d-shells and 
one in the O p-shell) on a Ni-O-Ni cluster. The dashed 
line signals a possible singlet between the O-hole with 
one of the Ni- holes producing a S=l/2 state as was de- 
scribed before. However, it is clear that the Ni partner of 
the oxygen hole can be easily switched from right to left 
producing an effective hopping process of the S=l/2 low 
energy state as shown in Fig. la. Other hoppings do not 
correspond to a mere interchange of the states but a spin 
flip can also occur. Exploring all possibilities, it can be 
easily shown that the allowed processes are those listed 
in Fig. lb (plus their spin reversed analogs). These hop- 
ping terms have been also derived by Zaanen and Oles 
p8| using more formal arguments in their complementary 
analysis of triplet holes moving in a S=l/2 background. 

The hopping processes (Fig. lb) can be written in a 
compact, explicitly rotationally invariant form. The full 
one-band Hamiltonian then becomes 

H = J^S i .S j -t^P ij (S i .S j + l/2). (2) 

<ij> <«) 

The first term is the Heisenberg interaction arising from 
Eq.(l) at half-filling and strong coupling. It only affects 
the S=l spins. Calculating J directly from Eq.(l) is dif- 
ficult and thus here we simply take the exchange from 
experiments i.e. J = 0.03eV. The second term in Eq.(2) 
acts only when spins 1 and 1/2 share a link. Py simply 
permutes the two spins. The notation S; = (S x ,S y ,S z ) 
is used in the hopping term to indicate that this spin 
operator can act over both S=l and 1/2 states de- 
pending on what spin is located at site i. Their ac- 
tion is the standard for a spin operator. For example, 

srs+K = !; m j z - -V 2 ) = \/2l m f = °; m j z = +V 2 )- 

Finally, it is natural to assume that the realistic cou- 
pling regime is t > J although we have not derived t and 
J from Eq.(l). @ 

An accurate analysis of one dimensional (ID) Hamil- 
tonians like Eq.(2) can be done using exact diagonaliza- 



tion (ED) techniques. |2(],^l[] With this method, we first 
calculated the dispersion of one hole, e(q), using model 
Eq.(2) (see Fig. 2a). e(q) is not symmetric with respect 
to 7r/2 contrary to systems with antiferromagnetic order. 
The minimum q m i n at t/J = 2.0 is at 7r/3 for the chain 
used in Fig. 2a. In the bulk limit we expect that q m in 
will move with continuity as a function of t/J. While 
the hole ground state at most q values has S=l/2, the 
states close to q = it have S=3/2 i.e. the band in Fig. 2a 
is actually made out of two energy levels with different 
spins. In the range 0.5 < t/J < 4.0, the bandwidth W is 
approximately 0.2 — 0.3J (Fig.2b)i.e. much smaller than 
the W ~ 2 J reported for holes in 2D S=l/2 antiferro- 
magnets. [|o| Several of these features could be observed 
in angle-resolved photoemission experiments. The pecu- 
liar dynamics of Eq. (2) is responsible for the abnormally 
small W. Such a large effective mass increases the ten- 
dency to hole localization, providing a natural explana- 
tion for the lack of true metallicity in Y2- x Ca x BaNi05 
and in doped 2D nickelates. 

Fig. 2c shows S(ir,ui) at a nominal hole density 
x = 1/12 - 0.083 and t/J = 2.0 calculated with ED. g] 
The qualitative agreement with the INS data reported 
for Y 2 - x Ca x BaNi0 5 [g at x - 0.04 - 0.10 is clear. The 
spectrum shows two main features: peak I caused by spin 
excitations in the S=l chain away from the hole, and 
peak II induced by the addition of a S=l/2 hole. To jus- 
tify these identifications, let us study the q-dependence of 
peaks I and II (Fig. 3a). Since S(q, oj) at finite x has low- 
intensity structure in addition to the two main peaks (see 
inset of Fig. 2c), the results shown in Fig. 3a are just rough 
estimations but we believe the main qualitative features 
are robust. Clearly peak I follows closely the position 
of the one magnon excitation observed in undoped S=l 
chains, while peak II is virtually dispersionless. Both 
peaks rapidly loose intensity when moving away from 7r, 
in agreement with INS results where only q = 7r pro- 
vided a large signal. |25| ] 

To better understand peak II let us study S(q, tu) for 
the undoped S=l chain with open boundary conditions 
(OBC) (Fig. 3b). At the ends of the chain, extended 
S=l/2 states are formed |24|] which arc coupled form- 
ing singlets (S) and triplets (T), becoming degenerate as 
N — > oo. Peak II in Fig. 3b corresponds to the transi- 
tion between the S and T states. We believe this transi- 
tion also contributes to S(7r,w) in Fig. 2c. We explicitly 
checked that the main contribution to peak II in Figs. 2c 
and 3b comes from spins located near the ends of the 
S=l chain. Note that on Fig. 3b it seems strange that the 
strength of peak II is relatively large since it should van- 
ish in the bulk recovering the Haldane gap ~ 0.4J (peak 
I). The reason is that the "size" of the S=l/2 end-states is 
~ 6 lattice spacings |24jj and thus only when N ^ 12 (or 
the density of hole defects is x <C 1/12) is that the contri- 
bution of the S-T transition will be negligible. This lead 
us to predict that for zinc doping y ~ 0.10, substantial 
weight inside the Haldane gap should be observed with 
neutron scattering as for the case of Ca doping (experi- 
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ments thus far have only been done for y = 0.04 pj). 

To understand the interaction between the S and T 
states with the hole consider first the unphysical limit 
tC J. In this case at low energies the dominant states 
are formed by the direct product of the S and T with 
the hole S=l/2. For a ring with N finite and even, the 
end-states are in a triplet state separated by a gap from 
the singlet combination. The triplet T couples to the 
S=l/2 hole forming states of S=l/2 and 3/2, denoted 
here as |1/2(T)) and 1 3/2) , respectively. These states are 
split by a finite t/J, and |1/2(T)) becomes the lowest en- 
ergy state for a dynamical hole. The singlet S forms a 
S=l/2 state with the hole (|1/2(S)) ) which is located at 
a higher energy. In Fig. 3c results are shown for |1/2(T)) 
and 1 3 / 2} at t/J = 0.1. If N is odd the positions of the 
quasi-degenerate sextuplet |1/2(T)) - |3/2), and the dou- 
blet |1/2(S)) are reversed. It is remarkable that increas- 
ing t/J to realistic values the energy difference between 
these states remain a fraction of J. For example, in Fig. 3c 
the energies of the |1/2(T)) and 1 3 / 2) states are also 
shown at t/J = 2.0. The low energy physics in S(q, oS) (at 
x > 1/12) is still dominated by the transitions between 
these two states. If x is reduced from ~ 0.08 the rela- 
tive importance of these states should rapidly decrease, 
but for the densities used in INS (x = 0.10) they are very 
important as already shown in Fig. 2c. 

Finally note that as t/J increases further, the kinetic 
energy becomes dominant and the hole prefers to be sur- 
rounded by a cloud of ferromagnetically aligned spins. 
On the N=12 chain, the hole state with S=3/2 becomes 
the ground state at t/J ~ 5.74. This coupling is not un- 
realistic and thus the high-spin states observed in the 
magnetic susceptibility || could be caused by the pres- 
ence of low energy S=3/2 Zhang- Rice holes. |^{| We also 
found that studying two holes on a 10 sites chain (den- 
sity x = 2/10 ~ 0.20), the remnants of the Haldane gap 
are difficult to identify in S(q, oS). This prediction can be 
tested experimentally increasing x in Y2- x Ca x BaNi05 
beyond the current limit x = 0.10. 

Summarizing, a model for mobile carriers in NiO 
chains was here proposed. S(q, u>) shows Haldane gap 
states in agreement with neutron scattering data. We 
believe that theories where holes are modeled as spin im- 
purities can be distinguished from the dynamical holes 
approach described here by measuring the optical con- 
ductivity. The present theory would predict a substantial 
weight inside the charge transfer gap upon doping since 
holes arc highly mobile and finite u> precursors of Drudc 
peaks should appear in severed NiO chains. 
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Figure Captions 

1. (a) Graphical representation of a hopping process 
as explained in the text; (b) Allowed hopping pro- 
cesses in model Eq.(2). The numbers in the ket 
represent the z-projection of the S=l and 1/2 spin 
states of Eq.(2) at two arbitrary sites i and j, re- 
spectively. Pijl/2 in Eq.(2) enforces the absence of 
the hopping |1, -1/2) -> | - 1/2, 1). The spin pro- 
jections only change by 1/2 or -1/2 at each site. 

2. (a) e(q) on a N = 12 chain with one hole and pe- 
riodic boundary conditions (PBC) using Eq.(2). 
The energy is measured with respect to the low- 
est energy state. Open (full) circles correspond to 
t/J = 0.5 (2.0): (b) Hole bandwidth W (energy dif- 
ference between the highest and lowest energy in 
e(q)) vs t/J on the N=12 chain; (c) S(ir, u) of model 
Eq.(2) with t/J = 2.0, N = 12, one hole and PBC. 
The meaning of peaks I and II is discussed in the 
text. In the inset S(7r/2,w) is shown for the same 
set of parameters. 

3. (a) Energy of the excitation with the highest inten- 
sity in S(q, u) (parameters as in Fig. 2c). The full 
squares denote peaks I and II as they evolve with q, 
with the error bars representing their width. The 
open squares are the one magnon peak position on a 
N = 14 undoped S=l Heisenberg model with PBC; 
(b) S(ir, u) for the undoped S=l Heisenberg model 
on a 12 sites chain with OBC. The meaning of the 
peaks is discussed in the text; (c) Energy vs q of 
the relevant low energy states (8 sites chain with 



PBC and one hole). Open squares and circles de- 
note the states (defined in the text) |1/2(T)) and 
1 3/2) , respectively, at t/J = 0.2. Full squares and 
circles correspond to |1/2(T)) and |3/2), respec- 
tively, at t/J = 2.0. 
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